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Abstract. The phase transition of the electroweak vacuum induced by a strong magnetic 
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Close to the lower critical magnetic field, a perturbative method developed by MacDowell 
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is shown that the triangular Abrikosov solution constitutes a local minimum of the energy 
provided M H > M z . 
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I Introduction 



In the cosmological scenario, the electroweak transition between the symmetric SU(2) ® 
U(1)y and the broken U(1) EM phases results from a temperature dependence in the coef- 
ficients of the Higgs potential. The transition occurs at a critical temperature T c and may 
involve the coexistence of phases, depending on the order of the transition. 

It has been shown by Ambj0rn and Olesen Jl], |2| that the phase transition can 
be induced at zero temperature by a large magnetic field. Their analysis was done with 
particular values of the coupling constants, for which the nonlinear field equations simplify, 
corresponding to Mjj = Mz where Mjj is the Higgs-boson mass. The transition from the 
broken phase with Higgs field <j) = (f>o to the symmetric phase with (f> = was found to 
take place gradually as the magnetic field increases from = My^/e to B c \ = M|/e.[] 
Between these values the field equations for A^, Zn, Wu, and 4> admit solutions with lattice 
periodicity in directions perpendicular to the magnetic field The solutions have been 
obtained numerically and studied in the range B C 2 < B < B c \ || . The emerging vacuum 
structure resembles that of the mixed state in a type-II superconductor (5| with the order 
parameter given by a density of W-boson pairs forming a zero-charge condensate. The 
distinctive feature is that the electroweak vacuum is paramagnetic, i.e. the magnetic field 
is enhanced by the W condensate, while a superconductor is diamagnetic. 

The phase transition resolves an old problem of the vacuum instability (6|, f/J of the 
Weinberg-Salam model for magnetic fields of mag nitude B c2 = M^/e ps 10 20 Tesla and 

beyond. The instability is a consequence of the large magnetic moment of the W bosons, 

1 With this choice of labels for the critical magnetic fields, the phase transition at B C 2 is qualitatively 
similar to the well studied transition at B C 2 in a type-II superconductor. In Refs. jlj [| the opposite labels 
were used. 
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which, when aligned with the magnetic field, can lower the energy sufficiently to make 
the W condensate energetically more favorable than the trivial vacuum. The W pairs also 
couple to and act as sources for nontrivial Z and Higgs fields. 

Explicit analytical solutions to the full, nonlinear problem are unknown, even for the 
case Mh = Mz considered in Refs. jj], |2|. Approximate solutions can be obtained, however, 
in certain limits. For B near the upper critical field B c \ (0 ~ 0) a perturbative solution 
was recently given by Olesen ||. In the vicinity of the lower critical field B C 2 (4> ~ ^o) a 
perturbative method was first suggested by Skalozub ||], who found that, in this regime, 
the solutions with lattice symmetry coincide with the Abrikosov solutions [|nj for a type-II 
superconductor near B C 2. In Skalozub's derivation, the interactions involving the Z and <j) 
fields were approximated with a local quartic interaction. 



In a recent paper [11], MacDowell and Tornkvist applied a perturbative method to 
the full Weinberg-Salam model for general values of the coupling constants and solved the 
field equations exactly to lowest order in (B — B C 2). In this approach, the interactions 
mediated by Z and (f) were accounted for in an effective, nonlocal quartic interaction. 
An investigation showed that, for Mg > Mz, the triangular^] Abrikosov lattice solution 
represents the energetically most favorable configuration. 

In this paper we shall derive analytical expressions for physical observables, such as 
the energy density, as a function of the geometry of the lattice solutions and of (B — B C 2). 
The geometry is specified by a complex parameter r. In the complex plane, r and the 
real number 1 span a fundamental parallelogram for an arbitrary two-dimensional simple 
lattice. 



2 In Ref. jll|] it was denoted "hexagonal" 
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The paper is organized as follows. In section |I|the perturbative method developed in 
Ref. [11] is reviewed. Physical quantities are shown to depend on the function |W| 2 where W 
is the perturbative solution for the W field. A general theorem of magnetic flux quantization 
is presented. For any solution of the nonlinear electroweak field equations where U(1)em 
gauge-invariant quantities possess a lattice symmetry in directions perpendicular to the 
magnetic field, it is shown that each lattice cell must enclose an integer number of quanta 
of abelian magnetic flux. 



Section III is devoted to mathematical properties of the Abrikosov solutions. A 
Fourier expansion of |W| 2 is derived and provides a fast-convergent representation with a 
simple dependence on the lattice parameter r. This representation is used in section ^ to 
derive explicit expressions for physical observables in terms of r. It is shown analytically 
that the triangular lattice solution constitutes a local minimum of the energy, provided 
M H > M z . 

II Electroweak Ginzburg-Landau Theory 



In the unitary gauge, the Weinberg-Salam lagrangean density jl2| leads to coupled equations 
for A^, Z^, and the real Higgs field (f> (see Ref. [| y] for details). The equations for 



and are 



D^F, 



cos 9 Zp, + sin 6 V - ig(W*W v - W$Wp)\ W» - \g 2 <?W v , (1) 



d»f^ = ie[d»{Wlw v -Wlw ll ) + WrtF tlv -W»Fl v ] , (2) 

where D» = d» + ig(A»smO + Z^cosO), F„ v = D^W V - D V W^ f„ v = d„A v - d v A^ 
Z^v = d^Z u — d u Z^ and e = gsinO. The field is subject to the constraint (TT 



D V W V = W V [^-Z V -d v \Tx<b 2 \ . (3) 
Vcosfc' ' 



We shall investigate static configurations of fields which include a magnetic field B 
in the z direction and where all fields are independent of the z coordinate. One can show 
that it is sufficient to consider components of the vector potentials in the x and y directions 
and that no inconsistent z or time components develop dynamically. The resulting problem 
is then two-dimensional. 

Define the spin projection states W± = (Wi =F iW2)/y/2. For B < B c2 = My^/e, 
Eqs. (|) and (g) have the trivial solution W M = 0, B = f\ = const., A= (B x r )/2, Z = 
and (f) = 4>q. When the magnetic field B exceeds B c2 , perturbations about zero in the 
linear combination W+ appear to become tachyonic in the field equation, Eq. (jl]), because 
of the spin interaction term igsmOf^W^. The trivial vacuum thus becomes unstable with 
respect to the production of pairs of W bosons with magnetic moments oriented along the 
magnetic field. Stability is restored by the cubic term in Eq. ([!]) and by the back reactions 
Z, <ft — <fto, B — B C 2, and A — A c2 for which W+ acts as a source. In Ref. |Tl| l it was shown 
that these back reactions are all of order |W+| . In particular, B in this order obeys the 
linear relation 



where H is a spatially uniform field. 

Using Eq. (|3|), one can also show that the suppressed component W_ is of order 



B-e\W+\ 2 = H 



(4) 



W+| 3 . Therefore, for an average magnetic field B above and near B c2 , we can write the 



equation for zero-energy eigenstates W+ as follows: 



(V - ieA) 




2 



2 



•2 



(5) 



where A is merely a notation for a vector potential such that V x A = Bz. 
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We recognize Eq. (0) as a generalization of the first Ginzburg-Landau equation of 
type-II superconductivity. That equation is the special case where 0{\W + \ Z ) is a positive 

2 

constant times and corresponds to a hamiltonian density with an effective local 

quartic interaction. In contrast, the effective hamiltonian of the electroweak problem is 



nonlocal. It was derived in Ref. [11] by expanding the exact hamiltonian up to second order 



2 

in |W+| including back reactions on the fields. The resulting effective hamiltonian is 

£{r ) = ^B 2 -(eB-M 2 r)\W+\ 2 + ~g 2 [sm 2 d \W + \ 4 + M^U(r )\W + \ 2 ] , (6) 

where 

U ^ ) = hJ d2r '[ K o(Mz\f-r'\)-K {M H \r-f'\)} \W+(r')f , (7) 

Kq is a Bessel function and Mz, Mjj are the masses of the Z and Higgs bosons. 

The interaction involving U(f ) can be interpreted as an effective long-range inter- 
action between W pairs, mediated by Z and Higgs bosons. For B > B c2 , Eq. (§) shows that 
the minimal energy occurs for a non-zero W+ field. In the case Mjj < Mz, one can find 
configurations which yield a negative quartic term. It is then necessary to go to higher or- 
ders of perturbation theory, which is beyond the scope of this paper. There are indications 
that an electroweak generalization of the Nielsen-Olesen vortex may be the stable solution 



in this mass regime [14]. This cannot be verified with a perturbative approach. 



If Mh > Mz, stability is ensured by the quartic interaction. Then there exist 
perturbative solutions of Eq. (||) for which [VF+I 2 is periodic on a fine-grain lattice of paral- 
lelograms and uniform on a macroscopic scale. They are known, from previous work in 



superconductivity, as Abrikosov flux lattice solutions [ 10 1 . 



In this paper, the solutions will be reconsidered in the context of the electroweak theory 
and physical quantities will be expressed as functions of the geometry of the lattice. This 
geometry is specified by the two lattice vectors ax and af or, equivalently, by the complex 
numbers a and or, a £ with the correspondence r = Re r x + Im r y. In the complex 
picture, |W+| 2 is referred to as doubly periodic with periods a and ar. 

Introduce the notation t r = Re r and 77 = Im r. For a given lattice, the cell side a 
and area A = a?Tj are dynamically determined by the value of the average magnetic field 
B through the following theorem, which is valid also for nonperturbative solutions. 

Theorem (Flux Quantization Condition) For field configurations where U(1) EM 
gauge-invariant quantities, such as fyi, Z v , (ft, \W-\ and W-/W+, are doubly peri- 

odic, the abelian magnetic flux that penetrates each lattice cell of area A is quantized and 
restricted to the values BA = 2nk/e , where the integer k is the common winding number 
of the phases of W+ and W- . 

Proof: Through integration by parts, Eq. (Q) can be written 
dPf^ = 2ied> J -(W} l W v -WlW il ) 

- ie [{rrwrfWv - wlip^Wp 

+ ie \(D u W^Wfj, — (DvWfj) . (8) 



JV 



Let W+ = |W+|e* x+ and W- = \ W-\e tx ~ . The last term on the right hand side then 



3 The theorem has previously been shown to hold for perturbative solutions |n]] and, nonperturbatively, 
for the special case of the electroweak theory where Mz = Mh jl5| • This proof was constructed in collabo- 
ration with MacDowell. 
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becomes 



A u + —d v ( X + + X-) 
2e 



jl°P = -2e\\W+\ z + \W. 

-2egco S e(\W + \ 2 + \W^ z )Z l/ -e(\W + \ z -\W.\ z )d u ( X+ -X-) • (9) 



The vector d u (x+ — X-) is invariant under lattice translation by virtue of the periodicity of 
W-/W+. Consequently, its line integral around the boundary of a parallelogram vanishes 
and, with the requirement that fields be single valued, 

fd v x + d X » = j>d v X-d X » = 2^ k integer. (10) 

Thus the phases of W + and W- have the same winding number k. Using Eq. (||) and, again, 
the periodicity of W—/W+, one can show that the middle term of Eq. (||) is an invariant 
vector under translation. When the expression for j* op , Eq. ©, is substituted into Eq. (|), 
the latter divided by 2e 2 (|W+| 2 + |W_| 2 ) can be rearranged in the form 

-A u = ±-d v ( X+ +x-)+jr , (ii) 

where j™ is an invariant vector under lattice translation. Therefore the integral § j™ v dx u 
around the boundary of a parallelogram vanishes, and we get 

flux = - j A v dx u = i- j d u { X+ + X-)dx v = ^ I (12) 



The periodicity condition on W-/W+ is required for the theorem to hold. It emerges 
naturally, if one assumes that the phases acquired under a lattice translation correspond to 
a gauge transformation, where the vector potential is form-invariant but expressed about 
the translated origin |l3[| . We remind ourselves that W- and W+ have the same U(1) EM 
charge and that W—/W+ therefore is a gauge invariant quantity. 
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From the flux quantization condition we find that the side a is determined by the 



relation 



2irk 
eBTj 



(13) 



It is convenient to redefine the problem in terms of coordinates, where the sides of a lattice 
cell have lengths 1 and |r|. The lattice is then specified by the sole parameter r, and we 
can impose t> > with no lack of generality. 



Define the dimensionless quantities B, p, V(p), W, e, and k by 



B 



eB 



, f=ap, V(p)=g 2 U(r), W + 



e, k 



kir 

Ti 



— ■ (14) 



With these units, the effective hamiltonian, Eq. (P), becomes 



e(p ) = \B 2 - {B - 1)|W| 2 + \ [ |W| 4 + V(p )\Wf 



(15) 



and the critical magnetic field is B C 2 = 1. From Eq. (Q) it follows that, in the new units, 
B - \W\ Z = B- |Wr = h is a uniform field. Substituting this relation into the effective 
hamiltonian we obtain the space-averaged energy density in terms of the average magnetic 
field B. 

1 _o _ 1 r 2 

(16) 



e= l -B 2 -{B-l)\W\ 2 + \ 



(|W| 2 ) +V{p)\W\ 2 



The lattice geometry does not fix the overall normalization of the solution W. For a given 
geometry, the normalization can be determined by minimizing the energy with a fixed 
average magnetic field B. The resulting condition is 



(B-i)iwr + 



(|Wh +V{p)\W[ 







(17) 
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Physical quantities can then be expressed in terms of B and the quantity 



^« , (IB. 

(inf) 



which is independent of the W normalization. One finds 



and 



, d») 

,_iB-_!2z£ , (20) 
2 2 1 + /? ' V ; 



/l = 1 + ( B - 1 ) 1 -^ • ( 21 ) 



Using Eqs. (0), (|TJ|), and (13) we can write 



B sin 



where 



and 



V{p,m) = -L JcPp' " K„(mp'J^) \W(p+ p' f (23) 



" W = TO <X = Z ' ff) ' P4) 
For Mjj > Afz we have j3 > and therefore 1 < h < B. The uniform field /i will 
remain frozen at ,6 C 2 = 1 if and only if Mh/Mz = 1, which is in agreement with the result 
of Ref. 1 15] for that mass ratio. 



This concludes the general description of electroweak Ginzburg-Landau theory. We 
have identified the physical quantities of interest. The purpose of this paper is to express 
them, analytically, as a function of the lattice geometry or, more precisely, of the lattice 
parameter r. In order to do so, we shall have to find a representation of the lattice solutions 



that will make possible an evaluation of the integral and averages in Eqs. (23) and (|2 
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Ill Expansions of Abrikosov Solutions 

The Abrikosov flux lattice solutions of Eq. (||) were first derived in the gauge A = Bxy. 
We prefer the cylindrically symmetric gauge A = {Bz x r ) /2, as it leads to a quicker and 
more elegant derivation. 

The lowest order perturbative solutions of Eq. (pf) are obtained by replacing 0(|W+| 2 ) 



with an effective mass term M 2 , (see Ref. [|TT[). After transforming to the rescaled units 
defined by Eq. fli~4|), and after expressing the vector p in cylindrical coordinates (p, ip), the 
equation becomes 

1 d d Id 2 „ . d o i 2k ( _ Mr 



W = . (25) 



It is easily checked that the angular momentum eigenstates with eigenvalues m : m > 

W m (p, <p) = -?L= exp(-i K p 2 ) p m e^ (26) 

Vvrm! 2 

and orthonormality condition 

d 2 p W m '(/9, if)* W m {p, if) = 5 m ' m (27) 



are infinitely degenerate solutions of Eq. (^5[) , satisfying the eigenvalue relation 

+ $ ■ (28) 
If we introduce the complex variable z = pe 1 ^, the most general solution is 

OO -i 

W(p, V) = Y1 c™W m (p, tp) = exp(--Kzz*) f(z) , (29) 

m=0 1 

where f{z) is an arbitrary analytic function. 
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i 1 2 

We are interested in solutions with |W| invariant under the lattice translations z — > z + 1 
and z — > z + r. Consider therefore the transformation properties of the theta function. 

0i(7r(z + l)|r) = -i?i(7rz|r) , 

i?i(7r(z + r)|r) = — exp(— iixr — 2iirz) -&i(tvz\t) . (30) 

The change of modulus under the second translation can be compensated for by attaching 
a prefactor with suitable transformation properties. If we make the choice 

f(z) = (2r 7 )3 exp(i K z 2 ) ^(ttz\t) (31) 

for / in Eq. (|29|), it is easily seen that W will transform by at most a phase under the 
two distinct translations, provided k = 7t/tj. According to Eq. (|l4|), this is the solution^ 
corresponding to a single quantum (k = 1) of magnetic flux per lattice cell: 

0iMt) . (32) 

The normalization is chosen so that, as we shall see, | W| 2 is coordinate covariant under the 
modular group and the spatial average | W| 2 is equal to one. The solutions with higher k 
are given simply by [W(p, <p)] . 

Since |W| 2 is doubly periodic, it can be expanded in a Fourier series in the coordinates 
(n, v) defined by z = u + vr. From Eq. ( |32"D and the series representation 

1 °° 2 

^ z \r) = - (-l)"g (n+ ^ e^ 2n+1)nz , (33) 

* n=— oo 

where g = e t7TT , |W| 2 can be written 

OO OO 2 2 

|W(«,U)| 2 = (2T 7 )2 V V (-l) n+ri ' e i7rr [ 1 '+(™+5)] g-WT^B-^'+i)] e i2ir(n-n>)u _ 
n=— oo n'=— oo 

(34) 

4 A proof of uniqueness of solutions was provided in Ref. . 
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W(p,<p) = (2r I )3exp 



7T 



— (z 2 
2r/ 



zz 



By trading the dummy index n for k = n — n 1 , the Fourier components in the u coordinate 
are already explicit. Since the expression is not termwise periodic in v, integration in v 
over merely a period will not help. Consider instead the continuous Fourier transform with 
respect to the v coordinate, 

\W(u;p)\ 2 = J dv e- i2npv \W(u,v)\ 2 

oo / 1 7 i2 \ oo 

= jT (-l) fc e^exp - T| " Pl \ e ^p(k+i) £ e i2* P n' _ (35) 

fc=— oo \ 1 I n'=— oo 



From the inverse Fourier transform and the Poisson formula 

oo oo 

E ei2npn '= E *(p-o . ( 36 ) 

n'=— oo /=— oo 

one then obtains the result 

oo oo / I t 7| 2 \ 

| W ( U , V )| 2 = £ E (-1) M+W exp ^ e «hr(fcu+to) (3?) 

k=— oo !=— oo \ / 

This representation converges extremely fast and uniformly on the plane. We remark 
that the expression is invariant under the modular group generated by the transformations 
r — > t' = r + 1 and r — > t' = — 1/r, and under reflexion in the imaginary axis r — > r' = — r*, 
provided u and u transform covariantly, i.e. n — > n', v — > where n + = n' + t)V [pTlf| . 
With the chosen normalization, the constant term 1 can be identified with the average value 



The Fourier expansion of |W| , Eq. (37), facilitates considerably the evaluation of 



integrals which occur in expressions for physical quantities in theories which allow Abrikosov 
flux lattice solutions, such as Type-II superconductivity or the electroweak theory. In 
particular it will enable us to write down an analytical expression for the parameter (5 [see 
Eq. (|l~8|) 1 in terms of the lattice parameter r. 
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IV Analytical Expressions for Physical Quantities 



In the lowest order of perturbation theory, physical quantities depend on the lattice geo- 
metry through dimensionless parameters which are specific to the Abrikosov solutions but 
independent of their overall normalization. In the theory of superconductivity, the naturally 
arising geometrical quantity is the Abrikosov parameter 



|W| 4 ~ / n\kr-lf \ 
^ = 2 = ^ ex P ~ • ( 38 ) 

(|W| 2 ) k,l=-oo \ Tl J 

The last equality was obtained by extracting the constant term in the Fourier expansion of 
| W | 4 that resulted from squaring Eq. (37). We remark that the right-hand side of Eq. (|3§1), 



and in fact any function of \kr — l\ /Yj, summed over all integers k and I, is modular 
invariant. 

According to section GIL the corresponding quantity in the electroweak theory is 



p _ v(p)\w\ 2 

The new feature here is the nonlocal quartic interaction in the numerator. It is described, 
as shown in Eqs. (^) and (p3[), by an integral kernel. 

In order to evaluate /3, we first find the Fourier expansion of the function V(p , to) 



defined by Eq. (23) and restrict ourselves to the case k = 1. With the representation 



„ , , 1 f°° dt ( x 2 \ 

Ko(l) = _y o _ exp _ ( -_ , (39) 



integrations are straightforward, and one finds 

oo 

V(u,v ;to) = (- 1 



k, l=—oo 



9 ir\kT — l\ 
to 2 + 2- 



Ti 



(40) 
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The average V(p ,m)\W\ is then obtained by multiplying Eqs. (40) and ( |37[ ) together 



and extracting the constant term. With mz and run defined by Eq. (p4j), the resulting 
expression for f3 is 

/3= KT,mz)-b(T,m H ) (41) 
B sin v 

where 



6(r,m) = 



E 



dwr) 

7r|A;T — l\ 



k, l=—oo 



m z + 2- 



Ti 



17 r 1 2 



exp 1 U . (42) 



The quantity (3 depends on, besides r, the masses of the two bosons that mediate 
the interaction and the redefined magnetic field B. It is therefore not a scale independent 
geometric parameter in the same sense as the Abrikosov number [3a- The dependence on B 
enters through the size of the flux lattice, which becomes significant with the introduction 
of interaction scales M^ 1 and Mjj 1 in the nonlocal kernel. We can, however, consistently 



set B = 1 in Eq. (41) within the order of perturbation theory we are considering. 

The above calculation can be done also for k > 1. Although the results are not 
as elegant and multiply convoluted sums abound, considerable computation time can be 
gained versus numerical integration. It has been shown numerically [|T|, |5| that the solutions 
with k > 1 have higher energy, and for this reason they are not the focus of this paper. 

The behavior of (3 in terms of the lattice parameter r has been investigated numer- 
ically in Ref. fll[| , where it was found that, for Mjj > Mz, the global minimum of (3 occurs 
at r = e* 71 / 3 corresponding to a triangular lattice. With the above results, it is now possible 
to show analytically that r = e l7r / 3 is a local minimum. 
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From Eqs. ( pl| ) and ( ]42| ) we have 

B sin 2 



1 fm H 

= _ . 2n / dm An , (43) 

m z 



where 



1 oo 

/3 m = -—b(T,m)= Yl E >n(M 2 ) , (44) 

A;, / = — OO 

fln(s) = 2 e-* , (45) 

(m/ + zx) 

and the rescaled lattice vectors 0J kl are given by uj m = ^J^(kr — I). The properties of j3 
can be demonstrated by analysis which holds true for each f3 m separately. Introducing 

d if d a \ a if d a \ 



one finds 



dr* 2 V 9t r 5tj / 3t 2\ dr R dr 1 



-Q^r = 1^ "hi F m (|w«| ) . (47) 

For the square (n = 4) and triangular (n = 6) lattices with n-fold rotational symmetry 
and r = e i27r / n , all non-zero lattice vectors u appear in n-tuplets {wr r , r = 1 . . . n}. Then, 
from the cancellation of phases within each n-tuplet, the right-hand side of Eq. (|47| ) sums to 
zero, and therefore the square and the triangular lattices are stationary points with respect 
to the variables t r and r z . This could have been shown directly from modular invariance 



In order to determine the type of local extremum that the triangular lattice consti- 
tutes, we must examine the eigenvalues of the Hessian matrix 

/ d 2 /3m d 2 (3 m \ ( d 2 /3 m d 2 /3m 

H= \ *** ^ UV ^ ^ \U , (48) 
\ 9^ drdT* ) \ dr R dri "dTT 5 

where U is a unitary matrix. 
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The elements of H are given by 

d 2 (3 m 



drdr* Atj 



— [ 2 M 2 F m '(\u k i\ 2 ) + \u kl \ A F m "{\uj M \ 2 )] , (49) 



k, l=—oo 



f) 2 B 1 °° r n 

= -472 E ^^^'(KH^+^^'C^I 2 )] , (50) 

k,l=— 00 

If the summation is performed separately over each n-tuplet {u;e 47rr//3 , r = 1 . . . 6} with 



common modulus, we see that the right-hand sides of Eqs. (^0|) and ( |5l| ) are zero by the 
cancellation of phases. The Hessian matrix is therefore diagonal with a double eigenvalue. 

To show that the eigenvalue is positive, examine the expression 2xF m '(x)+x 2 F m "(x) 
that occurs in each term of Eq. (|49|). It is easily shown to be positive for all m > 0, provided 
x > 2. Now for the triangular lattice we have 

hu| 2 > - = -7= > 2 (52) 
T i V o 



for each non-zero lattice vector t^. Therefore, the right-hand side of Eq. (49) is positive 
and, for all m > 0, it follows that j3 m has a local minimum at r = e* 71 ^ 3 . From Eq. ([43|) one 
then concludes that /3 has a local minimum for this value of r, provided Mh > Mz- 

The analysis of the Abrikosov parameter /3a can be carried out similarly. An inves- 
tigation of the perturbative properties of theories with more general quartic interactions is 
underway |lq|. 



According to Eq. (20), the energy is a monotonically increasing function of (3. If we 
assume that Mjj > Mz, the global minimum of occurs at r = e l7r / 3 Jll]], and it follows 
that the triangular lattice solution represents the ground state of the electroweak vacuum 
at magnetic fields above and close to the critical field 
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